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Abstract. A method is presented for the calculation of the one-body (IDM) and two-body (2DM) 
density matrices and their Fourier transforms in momentum space, that is consistent with the 
requirement for translational invariance (TI), in the case of a nucleus (a finite self-bound system). 
We restore TI by using the so-called fixed center-of-mass (CM) approximation for constructing an 
intrinsic nuclear ground state wavefunction (WF) by starting from a non-translationally invariant 
(nTI) WF and applying a projection prescription. We discuss results for the one-body (OBMD) 
and two-body (TBMD) momentum distributions of the '^He nucleus calculated with the Slater 
determinant of the harmonic oscillator (HO) orbitals, as the initial nTI WF. Effects of such an 
inclusion of CM correlations are found to be quite important in the momentum distributions. 


INTRODUCTORY REMARKS 


The last years the interest in the study of nuelei from both experimental and theoretieal 
point of view involves, besides the IDM 



( 1 ) 


and the OBMD 



( 2 ) 


also the 2DM 



( 3 ) 


and its Fourier transforms, for instanee the TBMD 




/ 


gip.(n-ri,)ei^-(f2-f2')p[2l(Pj^P2;rp,^20dVidVpd3r2dVy. (4) 


The above quantities provide among others information on the short-range eorrelations 
(SRC) in nuelei. 


A prominent role towards the experimental investigation of the 2DM and related 
quantities is played by the study of eleetromagnetieally indueed two-nucleon emission 
{Y,NN), {e^e'NN), carried out with high precision in photon facilities (ELSA, MAMI) 
and electron accelerators with high energy 100% duty-cycle beams (Jlab, MAMI) dlli. 

One of the theoretical issues still under discussion is a proper consideration of TI and 
the respective separation of spurious CM effects.The latter contaminate the calculated 
observables, when the independent-particle shell model and theories which take also 
dynamical correlations into account (e.g., Brueckner-Hartree-Fock, Variational Monte 
Carlo) are used, and inhibit the extraction of reliable information on the intrinsic proper¬ 
ties of nuclei directly from the experimental data (see, e.g., ilillll and Refs, therein). 

In this context, we prefer to deal with the intrinsic OBMD and TBMD and the as¬ 
sociated IDM and 2DM ||3]. They appear naturally when calculating the cross sections 
of nuclear emission in the plane-wave impulse approximation, i.e., neglecting the final- 
state interaction effects and the meson-exchange current contributions to the electro¬ 
magnetic interactions with nuclei. Eqs. (HJ)-® are modified by properly replacing the A 
space/momentum vectors by the respective Jacobi coordinates ivl]. In particular, we have 
calculated the expectation value for the TBMD 

^ A(A-l)('Ei„t|fiS(P:^)l'J'int), (5) 

where P = Ea=i Pa *^he total momentum operator for a nucleus of A nucleons and 'Lint 
is the intrinsic WE of the ground state. 


EVALUATION OF THE INTRINSIC OBMD AND TBMD 

The intrinsic quantities of interest are determined as the expectation values of appropri¬ 
ate A—particle operators, that depend on the Jacobi variables, in intrinsic nuclear states 
(see Eq. ®). An algebraic technique 0, [1] is applied for their evaluation, based on the 
Cartesian representation, in which the coordinate and momentum operators are linear 

combinations of the creation a' and destruction a operators for oscillating quanta in the 
three different space directions (see, e.g., 0). By following a normal ordering proce¬ 
dure, these intrinsic operators can be reduced to the form: an exponential of the set {a ^} 
times other exponential of the set {a}. For example, in the case of the intrinsic TBMD 
we get the representation 

T]g(p,k) = d3Aid3A2e-^^-^'e-‘^'^2£^^^(Xi,A2), (6) 

= e-^e-—^Di(z)...DA-2(?)OA-i(i*2)OA(i'i), (7) 

where da{y) = e“7''«ae“7«a (a = 1,... ,A), A = (Xi -f A2)/2, A = Xi —12 and 
-^1 = ^(^^-0) ’ ^2 = ^(^A+0) , z=-\^fA, 



with po the oscillator parameter in the momentum space. The Tassie-Barker-type factors 
in Eq. (|7]) appear in a model-independent way, i.e. they result from the intrinsic operator 
structure and do not depend on the intrinsic WF 'Tint, which is yet to be determined. 

In the following, we use the intrinsic unit-normalized WF constructed from a given 
nTI WF 'F, following Ernst-Shakin-Thaler (EST) (fixed-CM approximation) 

= (.R = 0|'F)/[('F|5(I)|'F)]1/2, (8) 

where |^) is an eigenvector of the CM operator R = and = 0) = 0. 

Here the bracket | ) is used to represent a vector in the space of the CM coordinates only. 


CALCULATION WITHIN THE INDEPENDENT PARTICLE 
MODEL: APPLICATION TO ^HE AND DISCUSSION 


One can show that if I'F) is the Slater determinant |Det) composed of single-particle 

(s.p.) orbitals ^,(a) (a = 1,... ,A), then the evaluation of the distribution is 

reduced to the evaluation of the matrix element 


(Det|di(z') ■ ■ ■ dA-2{zl)dA-i{x'2)dA{x'\)\T>a) = (Det(-v'i,-x'2,-z')|Det(x'i,x'2,z'))5 

that depends on some new complex vectors x'\,x' 2 ,z' 10] ■ The OBMD is evaluated in a 
similar way. The Slater determinant \T>e.i{x\,x’ 2 ,z’)) is deduced from the original |Det) 
via a substitution \^i{a)) Oal^iioc)). Further analytical evaluations are simplified 
with the HO orbitals. For example, in the simplest case of the configuration, which 
we encounter in the ^He nucleus, one can easily see that the matrix element is equal 
to unity since the s.p. state IO 5 ) coincides with the vacuum state of the Cartesian 
representation, viz., a|000) = 0, so that exp{a}|000) = |000). Moreover, the results 
obtained in this case are independent of the projection treatment used (see Eq. dHl)). 

In this way we have obtained the CM-corrected OBMD of "^He 

43/2^3 ^ 2,2 

VESlip) = VeSt(p) = ^3/2 ^3^2 ™ 

VS. the OBMD in the HO model without CM corrections 

»?<p(p) = >7.p(rt = (10) 


and the CM-corrected TBMD of "^He 



12 - 




^crrif 




-^k P'^^cm 


vs. the TBMD in the simple HO model ifioll 


( 11 ) 





The value for bQ {bcm = ■\/4/3Z?o) is obtained by equating the eharge rms radius of "^He, 
^rms = simple HO model + b^p in the HO model with 

CM corrections), with its experimental value (rrms =1-67 fm) and by taking the proton 
rms radius bp equal to 0.8 fm. We find bo, = 1.197 fm and bcm = 1.382 fm. 

From Eqs. (l^- (fT^ a shrinking of the distribution 7]est(/’) (^est(^’^)^ respect 

to risp(p) (T]|p (p,k)) follows, i.e., each of these distributions, after being CM-corrected, 
increases in its central but decreases in its peripheral region. More exactly, in the 
case of the TBMD this effect is related to the two-dimensional surface given by the 

function ri{p,k,cosy) = t 7 EST(^’^) variables p and k at each value of the angle 

7 between the vectors p and k. As shown in fl, the shrinking of the OBMD plays an 
essential role in getting a fair treatment of the data on the inclusive electron scattering 

in the GeV region. Another prominent feature of Tf^l^ip-ik) is its asymmetry due to 
the 7 —dependence. Fig. 1 demonstrates these changes for k = kpp, where kp is positive 
(negative) for k in the same (opposite) direction with respect to p. The H02 curves 
correspond to the calculations in Ref. m (bo replaced by bcm in Eq. (12)). Among 
the evident quantitative changes we observe the shift of the peak from kp = 0 towards 
negative k n’s, for p 7 ^ 0, due to a specific correlation induced by the CM fixation. 

In Ref. iUct the dimensionless quantity 

^(p,k) = p^^\p,k)/r](p)r](k) (13) 

was introduced as a measure of different correlations. In the complete absence of corre¬ 
lations ^ should be equal to 1 — 1 /A. For a finite, self-bound interacting fermion system, 
deviations of ^ from the above value is a measure of CM and(or) statistical and(or) dy¬ 
namical correlations. For the nucleus ^He in the simple HO (where A equals the level 
degeneracy of the only occupied state, thus statistical correlations are not active), 

^^1-1/A = 0.75 (14) 

for all p and k. After fixing the CM, ^ depends on p, k and 7 

&st(p,^) =0.89493e"5^"^-e-5^"^-e-i^-^*-. (15) 

In Fig. 2 logiQ^ is plotted as a function of cos 7 for selected values of p and k. ^ is 
significantly reduced in forward angles. The EST TBMD favors momenta of opposite 
directions as compared to the product of the two OBMD. The same holds for the TBMD 
of "^He if Jastrow-type SRC are included, as in Ref. (see the corresponding ^ for the 
case p = k = 4 fm, plotted in Eig. 2 (dotted line)). It is anticipated that within the EST 
approach additional corrections due to SRC will appear at high values of p and/or k and 
that they will be larger when p and k are antiparallel. 


PROSPECTS 

The method presented here is sufficiently flexible to be applied to a combined consider¬ 
ation of the CM and SRC. The latter, being introduced by means of Jastrow correlations. 
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FIGURE 1. The TBMD of for ^ ^ in the HO model with and without CM correlations (HOI; 
b - H02: b - bcra)- 
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FIGURE 2. The quantity logjQ ^ of “^He for the indicated values of p and A: as a function of cos 7 , where 
7 is the angle between p and k. Full, long and short-dashed lines: with CM correlations; dot-dashed: HO 
model without CM correlations; dotted: including SRC but not CM correlations (b = bcm)- 






do not violate the TI. We are also planning to extend our elaborations to other Z = N 
light-medium nuclei. This general formalism can be helpful in studying other two-body 
and many-body quantities. 

This research has been supported in part by the University of Athens under Grant No 
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